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Abstract 

Self-similar  multifractals  have  a  wavelet  transform  whose  maxima  define 
self-similar  curves  in  the  scale-space  plane.  We  introduce  an  algorithm  that 
recovers  the  affine  self- similarity  parameters  through  a  voting  procedure  in 
the  corresponding  parameter  space.  The  voting  approach  is  robust  with 
respect  to  renormalization  noises  and  can  recover  the  value  of  parameters 
having  random  fluctuations.  We  describe  numerical  applications  to  Can¬ 
tor  measures,  dyadique  multifractals  and  to  the  study  of  Diffusion  Limited 
Aggregates. 
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1  Introduction 


Some  important  classes  of  mnltifractals  are  composed  of  functions  which  are 
partly  self-similar.  The  applications  of  Iterated  Functional  Systems  have 
shown  the  importance  of  such  fractals  in  image  processing  [6].  Some  phys¬ 
ical  phenomena  such  as  fully  developed  turbulence  [12]  or  diffusion  limited 
aggregates  [1]  also  seem  to  exhibit  some  sort  of  self-similarity.  The  synthesis 
of  self-similar  fractals  is  a  relatively  simple  algorithmic  issue  but  the  analy¬ 
sis  of  the  affine  self-similarities  of  functions  or  measures  is  a  more  difficult 
inverse  problem.  For  physical  dynamical  systems  that  are  believed  to  have 
some  renormalization  properties,  it  is  particularly  important  to  have  pro¬ 
cessing  techniques  that  can  estimate  the  renormalization  parameters  from 
the  experiments.  We  introduce  a  voting  algorithm  that  can  recover  such 
renormalization  parameters  from  noisy  data. 

A  multifractal  might  have  different  renormalization  properties  at  differ¬ 
ent  scales.  The  necessity  to  decouple  the  multiscale  components  of  signals 
motivates  the  use  of  a  wavelet  transform  to  analyze  their  renormalization 
parameters.  The  renormalization  properties  of  a  function  appears  through 
self- similarities  of  its  wavelet  transform,  in  the  scale-space  plane.  For  many 
types  of  problems,  it  is  easier  to  analyze  the  properties  of  the  wavelet  trans¬ 
form  local  maxima  rather  than  the  values  of  the  transform  at  all  locations. 
In  particular,  we  proved  that  these  local  maxima  locate  and  characterize 
the  singularities  of  functions  [11],  The  wavelet  transform  maxima  define 
geometrical  curves  in  the  scale-space  plane.  The  self- similarities  of  a  func¬ 
tion  appear  as  geometrical  self- similarities  of  the  wavelet  transform  maxima 
curves.  We  thus  concentrate  on  these  local  maxima  to  analyze  the  renor¬ 
malization  properties  of  multifractals. 

The  multiplicative  properties  of  the  wavelet  transform  maxima  is  directly 
related  to  the  parameters  of  affine  self-similar  multifractals.  This  was  proved 
by  Bacry  et.  al.  [4]  in  the  case  of  Cantor  measures.  For  more  general  class  of 
multifractals,  they  conjectured  that  these  multiplicative  properties  is  related 
to  the  fractal  dimension  of  the  singular  support.  We  verify  this  result  for 
fractional  Brownian  motions.  Most  multifractals  encountered  in  physics 
or  image  processing  are  not  exactly  self-similar.  The  deviations  from  self¬ 
similarity  can  be  interpreted  as  a  renormalization  noise.  Voting  procedures 
such  as  the  Hough  transform  [5],  have  been  particularly  successful  to  recover 
parameterized  curves  from  noisy  data.  We  introduce  such  a  voting  scheme  to 
estimate  the  renormalization  parameters  of  the  wavelet  transform  maxima 
curves.  Section  4  gives  numerical  results  for  a  Cantor  measure,  a  dyadique 
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multifractal  function,  and  a  fractional  Brownian  motion.  An  application  to 
Diffusion  Limited  Aggregates  is  described  in  Section  5. 

2  Self-Similar  Multifractals 

An  interesting  class  of  multifractal  functions  is  solution  of  dilation  functional 
equations  of  the  type 


f(x )  =  J2pihf(h(x  ~  I'i)).  (1) 

2  =  1 

The  solution  might  not  be  a  function  but  a  distribution,  possibly  a  positive 
measure.  Throughout  the  article,  functions  and  measures  do  not  need  to 
be  differentiated  beyond  few  details  of  notations  that  are  left  to  the  reader. 
If  the  solution  f(x )  has  a  finite  non-zero  integral,  integrating  both  sides  of 
this  equation  yields 

n 

=  L 

2  =  1 

If  each  affine  renormalization  operation  is  locally  separable,  which  means 
that  for  1  <  i  <  n  there  exists  an  open  set  Dt  such  that  for  all  x  G  D; 

f{x)  =  Pilif{h{x  -  Ti)%  (2) 

the  function  is  locally  self-similar.  For  number  of  recurrence  equations,  f(x ) 
is  only  approximative^  equal  to  the  renormalized  version  on  the  right  hand- 
side  of  equation  (2),  which  introduces  what  we  call  a  renormalization  noise. 

Multifractal  Cantor  measures  satisfy  a  recurrence  equation  (1)  with  two 
separable  renormalization  components: 


f(x)  =  pihf{h(x  -  rq))  +  p2hf{h(x  ~  C >))■ 


The  renormalization  parameters  satisfy 


ri+2 17-'r2 


r^2f2 


1 
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Pl+  P2  =  1,  Pi  >0,  p2>  0. 

A  Cantor  measure  lias  a  support  included  in  the  interval  [— |,  |].  It  can  be 
constructed  recursively  as  follows.  We  define  a  uniform  measure  on  [— |,  |] 
and  divide  it  into  three  pieces:  a  uniform  measure  of  probability  pi  on  the 
interval  [—  |  a  zero  measure  on  the  interval  [—  |  |  —  ^],  and 

a  uniform  measure  of  probability  p2  on  the  interval  —  /2,  |].  The  left  and 
right  pieces  of  the  resulting  measure  are  then  subdivided  recursively  with  the 
same  renormalization  operations.  The  separability  property  (2)  is  satisfied 
for  the  domains  Di  =]  —  | |  +  ^-[  and  D2  =]|  —  |[. 

The  global  singular  properties  of  a  multifractal  is  characterized  by  the 
spectrum  of  singularity  /(a),  which  is  the  Hausdorff  dimension  of  the  set  of 
points  where  the  multifractal  has  singularities  of  strength  a.  For  a  Cantor, 
the  support  of  /(a)  is  an  interval  [tty,  0:2] 7  with  a.\  =  —  and  a2  =  — 

[9].  If  a  =  oq  =  a2,  then  all  singularities  have  the  same  strength  and  the 
Cantor  is  a  uniform  fractal  of  dimension  a. 

If  a  multifractal  is  self-similar  with  respect  to  affine  transformations  as 
in  equation  (2),  then  it  is  also  self-similar  with  respect  to  any  combination 
of  these  affine  transforms  and  their  inverse.  The  composition  of  n  affine 
transforms  whose  parameters  are  (pi,  1 1,  Tq)  with  m  transforms  of  parameters 
(p2 ihir2),  is  an  affine  transform  with  parameters  (p'lp'j' ■  Ci'/^S  rn_,n).  where 
rn,m  depends  upon  / 1,  l>-  cj  .  r2,  and  the  order  in  which  the  affine  transforms 
have  been  applied.  If  we  combine  as  well  the  inverse  of  these  transforms, 
we  obtain  the  same  result  for  n  £  Z  and  to  £  Z.  Our  goal  is  to  recover  the 
basic  affine  transforms  that  renormalize  any  given  multifractal  and  which 
yields  by  recombinations  all  other  affine  transforms  that  renormalize  this 
multifractal. 

3  Wavelet  Transform  Maxima  of  Multifractals 

The  wavelet  transform  is  a  natural  tool  to  analyze  the  multifractal  prop¬ 
erties  of  signals  since  it  separates  their  multiscale  components.  Arneodo, 
Bacry  and  Muzy  [2]  have  studied  its  application  to  several  types  of  physical 
multifractals.  A  wavelet  is  a  function  tp(x)  such  that 

f+00 

/  ip(x)dx  =  0.  (4) 
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The  wavelet  transform  of  a  function  at  the  abscissa  a  and  scale  s  is  defined 

by 

/H-OO  1  7*  qi 

f(u)-4{ - )du.  (5) 

-60  s  s 

The  wavelet  transform  of  a  measure  dp  is  defined  by  replacing  f(u)du  by 
dg(u).  Dilations  and  translations  of  functions  yield  simple  modifications  of 
their  wavelet  transform.  Let  g(x )  =  plf(l(x  —  a)),  we  derive  from  (5)  that 

Wg{s,  x )  =  plW f(ls,  l(x  -  r )).  (6) 

Let  us  suppose  that  our  wavelet  has  compact  support  included  in  the  interval 
\—K,  K].  If  a  function  satisfies  the  self-similarity  property 

Va-  e]a,  b[,  /(a)  =  plf(l(x  -  r)),  (7) 

then  for  any  s  such  that  5  <  ^4, 

Va  G  [a  +  Ks,  b  -  Ks ],  Wf(s,  a)  =  plW f(ls,  l(x  -  r)).  (8) 

The  self- similarity  of  /(a)  implies  a  self-similarity  of  its  wavelet  transform 
in  the  scale-space  plane.  However  the  converse  is  not  true.  If  the  wavelet 
has  n  vanishing  moments,  then  a  self-similarity  of  the  wavelet  transform 
implies  a  self- similarity  of  /(a)  up  to  a  polynomial  of  degree  n.  Indeed, 
for  such  wavelets,  the  wavelet  transform  ignores  polynomial  components  of 
degree  n.  Moreover,  the  wavelet  transform  might  be  self-similar  on  a  lim¬ 
ited  range  of  scales,  which  does  not  imply  a  self-similarity  of  /(a),  but  still 
proves  the  existence  of  a  partial  affine  self-similarity.  Hence,  the  wavelet 
transform  self-similarities  provide  a  finer  description  of  the  signal  renormal¬ 
ization  properties. 

The  wavelet  transform  is  a  two  dimensional  function  and  thus  represents 
a  large  amount  of  data  to  process.  To  characterize  the  singular  behavior  of 
functions,  it  is  sufficient  to  process  the  values  and  position  of  the  wavelet 
transform  modulus  maxima  [11],  These  modulus  maxima  have  been  used 
by  Bacry  et  al.  [4]  to  compute  the  spectrum  of  singularities  of  multifractals. 
Let  us  remind  that  a  wavelet  modulus  maxima  is  a  point  (sq,  x0 )  of  the  scale- 
space  plane,  where  |fT/(so,a:)|  is  locally  maximum  for  x  in  a  neighborhood 
of  a’0.  These  maxima  are  located  along  curves  in  the  scale-space  plane  (s,  a). 
Let  Mf(s,x )  be  the  maxima  map  defined  by  Mf(s,x )  =  Wf(s,x)  if  (s,a) 
is  a  modulus  maximum,  and  Mf(s,x )  =  0  otherwise.  Equation  (8)  implies 
that 

Va  G  [a  +  Ks,b—Ks],  M f(s,x)  =  plM f(ls,l(x  —  r)).  (9) 
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The  geometry  of  the  maxima  curves  are  thus  self-similar  in  the  scale-space 
plane.  Fig.  1(b)  shows  the  wavelet  transform  of  a  devil  staircase.  This 
function  is  the  primitive  of  a  Cantor  measure  built  with  parameters  pi  =  |, 
p2  =  /i  =  2  and  I2  =  4.  Fig.  1(c)  displays  the  position  of  the  modulus 
maxima  in  the  same  scale-space  plane.  The  self-similarities  of  the  wavelet 
transform  and  of  the  maxima  map  corresponds  to  the  self-similarities  of  the 
underlined  Cantor  measure.  The  wavelet  used  to  compute  this  transform 
is  shown  in  Fig.  2(a).  The  choice  of  wavelets  to  detect  and  characterize 
singularities  is  further  discussed  in  [11], 

As  a  consequence  of  equation  (9),  we  can  expect  that  the  number  of 
modulus  maxima  increases  when  the  scale  decreases.  Let  f(x )  be  a  Cantor 
measure  with  dilation  scales  1 1  and  I2.  Let  JV(s)  be  the  number  of  modulus 
maxima  of  the  wavelet  transform.  We  suppose  that  ip(x)  is  such  that  the 
number  of  maxima  is  bounded  at  any  scale: 


Vso  >  0,  3/i0,  Vs  >  s0,  JV(s)  <  K0.  (10) 


Bacry,  Muzy  and  Arneodo  proved  that  [4]  the  decay  of  N (s)  satisfies 

(11 

where  f3  is  given  by 

*£  +  /£  =  1-  (12) 

For  all  wavelets  used  in  practice,  the  constraint  on  the  finite  upper  bound 
(10)  is  always  satisfied.  Equation  (11)  proves  that  for  Cantor  measures  the 
decay  of  the  number  of  maxima  depends  upon  the  multiplicative  properties 
of  the  measure  characterized  by  the  dilation  scales  /  [  and  /2.  When  the  Can¬ 
tor  measure  is  uniform,  which  means  that  all  singularities  have  a  strength 
equal  to  a  =  — since  pi  +p2  =  1,  we  can  derive  that  / 3  =  —a. 
For  uniform  Cantors,  the  decay  rate  of  the  number  of  maxima  is  thus  equal 
to  their  fractal  dimension.  For  general  Cantors,  —f3  is  the  fractal  dimension 
of  the  support  of  the  measure. 

Fig.  3  gives  the  number  of  maxima  across  scales  corresponding  to  the 
wavelet  maxima  map  shown  in  Fig.  1(c).  If  we  ignore  constant  factors,  the 
devil  staircase  has  the  same  renormalization  properties  as  the  underlined 
Cantor  measure.  A  linear  regression  on  the  number  of  maxima  in  Fig.  1(c) 
yields  3  =  —0.70  and  indeed 

2P  =  0.99. 


lim  = 

log  S 
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Figure  1:  (a)  Primitive  of  a  Cantor  measure  constructed  with  pi  =  2/3, 
p2  =  1/3,  l\  =  2,  and  1 3  =  4.  (b)  Wavelet  transform  of  the  Devil  Staircase 
shown  in  (a),  (c)  Modulus  maxima  of  the  wavelet  transform  shown  in  (b). 


(a)  (b) 

Figure  2:  (a)  Graph  of  a  wavelet  ip(x)  with  compact  support  and  one  van¬ 
ishing  moment.  It  is  a  quadratic  spline,  (b)  Graph  of  the  Primitive  0(x ) 
with  compact  support.  It  is  a  cubic  spline. 


Figure  3:  This  curve  gives  the  decay  of  loglV(s)  as  a  function  of  logs,  where 
N(s)  is  the  number  of  wavelet  transform  modulus  maxima  at  the  scale  s  of 
the  Devil  Staircase  shown  in  Fig.  1(a). 
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For  more  general  multifractals,  Bacry  et.  al.  [4]  conjectured  that  the 
number  of  maxima  N(s)  decays  like  s~d  where  d  is  the  fractal  dimension 
of  the  singular  support  of  the  multifractal.  The  next  theorem  verifies  this 
results  for  fractional  Brownian  motions.  Fractional  Brownian  motions  are 
Gaussian  processes  that  are  self-similar  in  a  stochastic  sense,  with  a  power 
spectrum  proportional  to  \^nr+ r )  for  some  exponent  1  >  H  >  0.  Their 
realizations  are  not  affine  self-similar  functions  but  fractal  curves. 

Theorem  1  The  expected  density  D(s)  of  modulus  maxima  of  the  wavelet 
transform  of  a  fractional  Brownian  motion  satisfies 

D(s)  =  -, 
s 

where  C  is  a  constant  that  depends  only  upon  the  wavelet  tp(x)  and  the 
fractional  exponent  H . 

The  proof  of  the  theorem  and  the  value  of  the  constant  C  are  given  in 
the  Appendix.  The  decay  of  the  number  of  maxima  for  fractional  Brownian 
motions  is  thus  inversely  proportional  to  the  scale.  Since  fractional  Brown¬ 
ian  motions  are  almost  everywhere  singular,  the  fractal  dimension  of  their 
singular  support  is  d  =  1,  so  this  theorem  verifies  the  conjecture  of  Bacry 
et.  al.  [4], 

4  Vote  for  Renormalization  Parameters 

Since  the  self-similarities  of  a  function  implies  geometrical  self- similarities  of 
the  wavelet  transform  maxima,  we  concentrate  on  these  maxima  to  find  the 
values  of  the  renormalization  parameters.  We  can  analyze  the  geometrical 
properties  of  the  maxima  lines  in  the  scale-space  plane.  This  approach  was 
also  studied  independently  by  Arneodo,  Bacry  and  Muzy  [3],  who  devel¬ 
oped  an  algorithm  based  on  the  bifurcation  points  of  the  wavelet  transform 
maxima,  in  order  to  detect  and  characterize  renormalization  maps.  In  most 
cases,  multifractals  are  approximatively  self-similar,  which  means  that  we 
must  take  into  account  renormalization  errors.  Many  algorithms  that  esti¬ 
mate  the  parameters  of  noisy  curves  in  an  image  plane  are  voting  procedures 
in  the  corresponding  parameter  space  [5].  We  describe  such  a  voting  algo¬ 
rithm  that  can  recovers  non-exact  renormalization  properties  in  multifractals 
by  using  the  information  provided  by  all  the  wavelet  transform  maxima. 


Affine  renormalizations  are  characterized  by  the  three  parameters  ( p ,  /,  r ) 
of  equation  (7).  We  introduce  a  voting  scheme  in  the  three  dimensional  pa¬ 
rameter  space  ( p,l,r ).  Renormalization  parameters  are  then  identified  as 
points  of  high  votes  in  the  space  (p,  /,  r ).  As  mentioned  in  the  previous  sec¬ 
tion,  if  the  multifractal  is  invariant  under  two  affine  transformations  specified 
by  (pi,h,ri)  and  (p2,  h,  ^2),  it  is  also  invariant  with  respect  to  a  family  of 
affine  transforms  specified  by  (pjp™,  rnim),  where  n  and  m  are  inte¬ 
gers.  This  produces  spurious  high  vote  peaks  at  the  corresponding  locations 
of  the  parameter  space.  To  obtain  a  regular  distribution  of  these  peaks,  we 
use  a  logarithmic  scale  for  p  and  /.  In  the  remaining  of  the  paper,  log  a-  is  a 
logarithm  base  2  of  a.  The  parameter  space  is  thus  represented  by  a  three 
dimensional  array  indexed  by  (logp,  log  /,  r),  called  accumulator  array.  Each 
bin  of  this  array  corresponds  to  a  cube  in  the  space  (logp,  log  I,  r)  of  size 
A  logp  X  A  log  I  X  Ar.  The  values  of  A  logp,  A  log!  and  Ar  depend  upon 
the  desired  precision  when  measuring  the  renormalization  parameters. 

Equation  (6)  proves  that  a  renormalization  by  (p,  1,  r)  maps  any  wavelet 
maxima  located  at  (si,  aq)  to  a  maxima  located  at  (s2  =  Isi,  x2  =  l(x  1  —  r)). 
If  /( x)  is  invariant  with  respect  to  the  affine  transform  of  parameters  (p,  /,  r), 
equation  (8)  proves  that 

P  lWf(s2,x  2) 

Any  maxima  in  the  neighborhood  of  (sj,  xi)  is  also  mapped  to  a  maxima 
in  the  neighborhood  of  (s2tX2)  with  the  same  affine  transformation.  Let 
(si,  xi  +  Ai)  be  the  location  of  the  closest  maxima  to  (si,  xi),  with  Ai  >  0. 
This  maxima  is  mapped  to  (s2,  x2  +  A2)  with 

A  2  _ 

Ai  s1 

For  any  pair  of  wavelet  transform  maxima  (si,  xi)  and  (s2,  x2),  if  the  closest 
left  maxima  of  (gj,  X|)  and  (.s2,  x2),  located  respectively  at  (si,  xq  +  Aj)  and 
(s2,  x2  +  A2)  satisfy  the  constraint  (14),  we  then  vote  for  the  renormalization 
parameters 


log  1  =  log!  — ), 

Sl 

(15) 

i  i  ,  Wf{si,xi) 

(16) 

Sl 

r  =  x  1  -  x2 — . 

^2 

(17) 

(13) 
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The  voting  algorithm  proceeds  as  follow.  The  bins  of  the  parameter  space 
array  are  first  initialized  to  0.  For  each  pair  of  wavelet  transform  maxima 
(si,  xi)  and  (s2,  x2)-  with  si  <  s2,  if  their  closest  right  maxima  satisfy  equa¬ 
tion  (14),  we  add  1  to  the  bin  (log  p,  log  f,  r)  defined  by  equations  (15-17). 
This  algorithm  yields  high  votes  in  the  bins  corresponding  to  the  renormal¬ 
ization  parameters  of  the  multifractal.  The  renormalization  noise  can  spread 
the  votes  for  a  given  set  of  renormalization  parameters  across  several  bins  of 
the  accumulator  array.  The  bin  size  Alogp,  A  log  I,  Ar  of  the  accumulator 
array  must  therefore  be  adapted  to  the  amount  of  renormalization  noise  in 
order  to  avoid  splitting  votes.  Once  the  vote  is  done,  we  select  the  parame¬ 
ter  indexes  with  high  votes,  where  the  value  of  the  vote  is  locally  maximum 
in  a  three-dimensional  neighborhood  of  the  accumulator  array.  If  the  algo¬ 
rithm  is  successful,  the  peaks  of  highest  votes  provide  the  renormalization 
parameters  of  the  multifractal.  Instructions  to  obtain  a  copy  of  the  software 
implementing  this  algorithm  are  available  through  anonymous  ftp  at  the  ad¬ 
dress  cs.nyu.edu,  in  the  file  README  of  the  directory  /pub/wave/software. 

We  first  illustrate  the  results  of  this  algorithm  for  the  devil  staircase 
with  1-l  —  2,  l2  =  4,  p1  —  |,  p2  —  |  and  r1  =  r2  =  f,  shown  in  Fig. 
1(a).  The  voting  is  based  on  the  wavelet  transform  maxima  displayed  in 
Fig.  1(c).  The  highest  votes  in  the  three  dimensional  accumulator  array  are 
at  (log pi ,  log /j,  rq)  and  (logp2,  log  /2,  r2).  Since  a  Cantor  set  is  also  invari¬ 
ant  by  affine  transforms  of  parameters  (PiP™,  rnm)  for  ( n,m )  G  Z2, 
the  algorithm  also  yields  peaks  of  votes  at  (nlog(pi)  +  mlog(p2),  nlog(/i)  + 
mlog(/2), rnm).  In  the  plane  (logp,logI)  these  votes  are  located  on  a  uni¬ 
form  grid  whose  interval  is  specified  by  (logpi, log li)  and  (logp2, log /2).  It 
allows  us  to  verify  that  the  multifractal  only  has  two  sets  of  renormalization 
parameters  whose  values  are  specified  by  the  locations  of  the  two  highest 
votes  in  the  parameter  space.  Fig.  4(a)  and  Fig.  4(b)  show  the  votes  in 
the  planes  of  coordinates  r  =  —0.26  and  r  =  0.38,  in  the  three  dimensional 
parameter  space  (logp.logf,  r ).  The  two  high  peaks  are  the  highest  votes  in 
the  three  dimensional  parameter  space. 

Dyadique  renormalization  equations  also  define  interesting  multifractals 
studied  by  Deslauriers  and  Dubuc  [7],  as  well  as  Daubechies  and  Lagarias 
[10].  The  function  shown  in  Fig.  5  is  the  solution  of  the  functional  equation 

/(x)  =  |/( 2(x  -  I))  +  /( 2x)  +  |/( 2(x  +  1)1,  (18) 

analyzed  by  Daubechies.  When  applying  the  voting  procedure  to  the  max¬ 
ima  of  the  wavelet  transform  of  f(x ),  surprisingly  there  are  4  peaks  that 
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Figure  4:  Votes  based  on  the  wavelet  modulus  maxima  for  the  Devil  Stair¬ 
case  of  Fig.  1(b).  These  figures  give  the  number  of  votes  in  the  two  planes 
of  the  parameter  space  (logp,  log  /,  r),  where  the  highest  number  of  votes  are 
concentrated.  The  plane  (a)  corresponds  to  r  =  —0.26  and  the  plane  (b)  to 
r  =  0.38.  The  location  of  the  two  high  peaks  provide  the  renormalization 
parameters  of  the  underlined  Cantor  measure. 
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t'(x) 


Figure  5:  Solution  of  the  functional  equation  / (x)  =  l/3/(2(x  —  1/2))  + 
/(2s)  +  2/3/(2(as  +  1/2)). 

emerge  strongly  from  the  accumulator  array,  for  parameters  (p, ,  /,  r)  respec¬ 
tively  equal  to  (|,  2,0),  (|,  2,0),  (|,2,|)  and  (|,2,—  |).  Fig.  6(a-c)  show 
the  votes  in  the  planes  corresponding  respectively  to  r  =  —  |,  r  =  0  and 
r  =  |.  Beyond  the  four  peaks,  we  can  see  other  local  peaks  of  smaller 
amplitude  that  correspond  to  recombinations  of  the  four  main  affine  renor¬ 
malizations.  It  is  surprising  to  see  four  different  affine  self-similarities  since 
the  recurrence  equation  (18)  has  only  three  terms  and  the  parameters  can 
not  be  derived  directly  from  the  coefficients  of  the  equation.  However,  a 
more  careful  analysis  of  the  solution  f(x )  reveals  that  these  affine  trans¬ 
forms  are  indeed  the  true  renormalization  parameters.  Let  us  decompose 
f(x)  into  f(x)  =  g(x)  +  h(x),  where  the  support  of  h(x)  and  g(x)  are  re¬ 
spectively  [—1,0]  and  [0, 1].  One  can  prove  that  h(x )  is  the  primitive  of  a 
positive  Cantor  measure  centered  at  —  1  with  parameters  If  =  2,  /2  =  2, 
pi  =  1/3  and  p2  =  2/3,  and  g(x)  is  the  primitive  of  a  negative  Cantor  mea¬ 
sure  centered  at  |  with  parameters  1 1  =  2,  /2  =  2,  pi  =  1/3  and  p2  =  2/3. 
In  this  case,  the  voting  algorithm  did  indeed  reveal  non-trivial  self-similar 
structures  of  the  function  f(x). 

As  a  last  example,  we  describe  the  results  of  this  voting  algorithm,  when 
applied  to  a  fractional  Brownian  motion.  The  realizations  of  these  stochastic 
multifractals  do  not  have  affine  self-similarities,  and  indeed  do  not  yield  any 
dominating  peak  in  the  parameter  space.  To  visualize  better  the  distribution 
of  votes  in  the  parameter  space,  for  any  given  r  we  sum  the  votes  for  all 
values  of  p  and  /.  The  resulting  number  of  votes  as  a  function  of  r  is  given 
in  Fig.  7(a).  As  expected,  this  distribution  is  uniform,  which  indicates  that 


12 


vote 


vote 


o  -3 


(b) 


vote 


(c) 

Figure  6:  Votes  based  on  the  wavelet  transform  modulus  maxima  of  the 
dyadique  multifractal  shown  in  Fig.  5.  These  figures  give  the  number  of 
votes  in  the  three  planes  of  the  parameter  space  (logp,  log /,  r),  where  the 
highest  number  of  votes  are  concentrated.  The  planes  (a),(b)  and  (c)  cor¬ 
respond  respectively  to  r  —  —  r  —  0  and  r  =  The  location  of  the  four 
high  peaks  provide  the  renormalization  parameters  of  this  multifractal. 
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Figure  7:  Distribution  of  votes  for  a  fractional  Brownian  motion  with 
H  =  0.2.  (a)  Repartition  of  votes  as  a  function  of  r,  after  summing  over  all 
values  of  log  I  and  log  p.  (b)  Rep  art  ion  of  the  votes  in  the  plane  (logp,logZ), 
after  summing  over  all  values  of  r.  Brighter  pixels  correspond  to  larger 
votes.  There  is  no  high  peak  in  (a)  and  (b). 
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there  is  no  remarkable  translation  parameter.  For  any  given  (logp,log/), 
we  also  sum  the  votes  for  all  r  values  in  the  accumulator  array  and  obtain 
a  distribution  in  the  (logp,logI)  plane,  shown  in  Fig.  7(b).  Higher  votes 
correspond  to  brighter  points  in  the  image  plane.  No  particular  peak  emerges 
from  this  plane,  but  the  vote  has  a  clear  distribution  around  a  straight  line 
whose  equation  is  logp  =  /3 log /,  with  f3  =  —1.2.  The  slope  f3  is  related 
to  the  fractal  parameter  H  =  0.2  of  this  fractional  Brownian  motion,  by 
/3  =  —  1  —  H .  This  can  be  justified  with  the  following  intuitive  argument. 
The  algorithm  votes  for  a  parameters  p  in  equation  (13)  which  is  equal  to 
the  ratio  of  the  wavelet  transform  at  two  maxima  locations,  divided  by  l. 
Flandrin  [8]  proved  that  at  any  point  (si,*i),  the  variance  of  the  wavelet 
transform  of  a  fractional  Brownian  motion  satisfies 

E{\WB{s1,x1)\2)  =  Cs12H, 

where  C  is  a  constant  independent  of  (aq,  s L ) .  Hence  the  ratio  of  the  wavelet 
transform  at  two  maxima  locations  at  the  scales  iq  and  S2  is  on  average  close 
to  ( =  l~H ■  The  value  of  logp  will  thus  on  average  be  of  the  order  of 
{—H  —  l)log/.  This  explains  the  higher  density  of  votes  around  the  line 
of  equation  logp  =  (—  1  —  if)  log  I.  In  this  example,  the  voting  algorithm 
indicates  that  this  fractal  has  no  affine  renormalization  but  we  can  still 
measure  the  fractional  exponent  H  from  the  repartition  of  the  votes  in  the 
( log p, log!)  plane. 

5  Diffusion  Limited  Aggregates 

The  analysis  of  the  properties  of  fractal  growth  phenomena  remains  mostly 
an  open  problem  in  physics.  Diffusion  Limited  Aggregates  can  be  modeled  as 
a  growth  obtained  by  successive  accretion  of  random  walker  at  the  periphery 
of  the  cluster.  Arneodo  et.  al.  [1]  have  recently  introduced  techniques  based 
on  the  wavelet  transform  modulus  maxima  to  analyze  the  fractal  properties 
of  these  DLA.  Fig.  8(a)  shows  the  center  region  of  a  DLA  cluster  that 
contains  a  total  of  108  particles.  To  analyze  the  properties  of  such  clusters, 
Arneodo  et.  al.  [1]  have  constructed  a  measure  that  indicates  the  location 
of  intersection  points  of  such  a  cluster  with  a  circle  of  fixed  radius.  At 
intersection  points,  the  measure  is  equal  to  1  and  0  at  other  locations.  Fig. 
8(b)  such  a  measure  built  at  the  intersection  of  a  circle  and  a  DLA  cluster. 

The  wavelet  transform  is  normally  defined  with  respect  to  a  wavelet 
function  whose  integral  is  zero.  However,  to  analyze  positive  measures  of 
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compact  support,  one  can  use  a  wavelet  primitive  whose  integral  is  non-zero, 
and  obtain  the  same  properties  as  a  standard  wavelet  transform.  Let  4,1(x) 
be  the  primitive  of  compact  support  of  a  wavelet  ip{x)  °f  compact  support. 
The  function  /01(x)  might  have  a  non-zero  integral  and  thus  might  not  be 
a  wavelet.  Let  f(x)  be  a  positive  measure  of  compact  support  and  g(x )  be 
its  primitive  of  compact  support.  Let  us  denote  W1  f(s,x)  the  “pseudo” 
wavelet  transform  of  f{x)  computed  with  ip1{x)  by  using  the  same  formula 
(5)  as  a  wavelet  transform.  Let  Wg(s,x)  be  the  wavelet  transform  of  g(x ) 
computed  with  the  wavelet  ip(x).  An  integration  by  part  of  equation  (5) 
yields 

W1f(s,  x)  =  -Wg{s,x). 
s 

Hence,  the  pseudo  wavelet  transform  of  / (x)  computed  with  -01(a;)  has  the 
same  properties  as  the  wavelet  transform  of  its  primitive  computed  with 
respect  to  ip(x).  To  analyze  the  properties  of  DLA  measures,  we  use  the 
wavelet  primitive  tp1( x )  shown  in  Fig.  2(b).  Arneodo  et.  al.  [1]  showed 
that  the  number  of  maxima  of  DLA  measures  at  scales  2.2”  is  a  Fibonacci 
sequence.  One  interpretation  is  that  each  branch  of  the  DLA  splits  into  two 
branches,  one  larger  than  the  other,  which  then  splits  recursively  with  the 
same  procedure.  It  is  however  not  clear  how  to  relate  precisely  results  on  one¬ 
dimensional  DLA  measures  to  the  properties  of  two-dimensional  aggregates. 
In  this  section,  we  analyze  the  properties  of  DLA  measures  from  the  results 
of  the  voting  algorithm  described  in  section  4.  The  experiments  are  done  on 
50  DLA  measures  provided  by  Arneodo.  Fig.  9  gives  the  number  of  wavelet 
modulus  maxima  jV(s),  computed  with  50  DLA  measures,  for  .s  £  [22 , 28] . 
A  linear  regression  yields 

s 

log  s 

To  relate  this  result  to  Fibonnacci  sequences,  let  us  recall  that  a  Fibonacci 
sequence  (un)ne jq  satisfies 

linr  lQg^re-  =  log(»,  (20) 

n— >+oo  n 

where  <f>  =  1+2v''T  Equation  ( 19)  and  (20)  implies  that  the  number  of  maxima 
N{s)  may  be  a  Fibonacci  sequence  along  the  scales  s  =  a”,  if  a  =  which 
yields  a  =  2.13.  This  result  is  consistent  with  the  fact  that  Arneodo  et.  al. 
found  that  the  number  of  maxima  at  scales  2.2"  vary  almost  like  Fibonacci 
sequences. 


=  -0.63  ±  0.01. 


(19) 
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Let  us  now  study  the  renormalization  properties  of  these  DLA  with  the 
voting  algorithm.  To  simplify  the  analysis  we  shall  concentrate  on  the  pa¬ 
rameters  p  and  /  and  eliminate  the  translation  parameter  r.  We  make  a 
vote  for  the  parameters  ( log p,  log!)  from  the  wavelet  transform  maxima 
computed  with  the  50  DLA  measures.  Fig.  10(a)  shows  the  distribution  of 
the  highest  votes  in  the  plane  (logp,  log/).  The  votes  are  distributed  along 
a  line  of  equation  logp  =  —a  log/,  with 

a  =  0.63  ±  0.01. 

This  indicates  that  the  singularities  of  strength  a  =  0.63  strongly  dominate 
all  other  singularities.  It  is  interesting  to  observe  that  this  exponent  is  the 
same  as  the  one  measured  in  equation  (19)  from  the  decay  of  the  number 
of  maxima.  We  saw  in  section  3  that  such  a  situation  occurs  for  uniform 
Cantor  measures.  For  example,  a  uniform  Cantor  with  parameters  l-y  =  2.13, 
Z2  =  (2.13)2,  pi  =  0.618  and  P2  =  0.382  does  satisfy  this  property  and  the 
number  of  maxima  at  the  scales  an  for  a  =  =  2.13  is  a  Fibonacci  sequence. 

This  Fibonacci  Cantor  is  built  by  recursively  dividing  a  uniform  measure  into 
two  uniform  measures.  One  has  a  support  times  smaller  and  the  other 
one  a  support  / 2  =  l\  times  smaller.  When  the  branches  of  a  DLA  splits, 
Arneodo  observed  that  the  larger  branch  has  one  chance  out  of  2  to  be  on 
the  right  or  on  the  left  of  the  smaller  branch.  To  take  into  account  this 
fact,  one  can  build  a  simple  stochastic  Fibonacci  Cantor  model  where  the 
sub-component  of  size  is  located  randomly  on  the  left  or  on  the  right 
of  the  component  of  size  /2,  when  recursively  building  the  measure.  Fig. 
10(b)  displays  in  the  (logp,log/)  plane  the  distribution  of  votes  computed 
from  the  wavelet  transform  maxima  of  several  realizations  of  this  stochastic 
Fibonacci  Cantor.  As  expected,  the  votes  are  distributed  around  a  line  of 
equation  logp  =  —  aTog/,  with  a  =  0.63. 

For  any  fixed  log/,  Fig.  11(a)  and  11(b)  show  the  distributions  ofvotes  as 
a  function  of  log  /  when  summing  the  votes  along  all  logp  for  DLA  measures 
and  for  stochastic  Fibonacci  Cantors.  The  vote  distribution  of  stochastic  Fi¬ 
bonacci  Cantor  has  three  clear  peaks  at  log  /  =  1.1,  log  /  =  2.2  and  log  /  =  3.3 
which  are  multiples  of  the  basic  renormalization  parameters  log/i  =  1.1  and 
log/2  =  2.2.  For  the  DLA  measures,  there  are  two  peaks  at  log  /  =  1.0  and 
log  /  =  2.1  which  are  much  weaker  and  a  third  rebound  at  log  /  =  3.3,  which 
indicates  that  this  simple  stochastic  Fibonacci  Cantor  model  do  not  reflect 
precisely  all  structures  of  the  DLA.  One  can  refine  the  Cantor  model  by  let¬ 
ting  the  left  and  right  components  of  the  measure  have  positions  that  vary 
randomly  when  building  it  recursively,  but  it  is  also  likely  that  the  multi- 


17 


plication  process  involved  in  DLA  has  a  more  complicated  structure  than  a 
Cantor  subdivision.  A  more  precise  analysis  of  the  distribution  of  the  votes 
of  DLA  measures  that  includes  the  translation  component  r.  can  help  specify 
this  initial  Fibonacci  model,  but  a  detailed  analysis  of  the  physics  of  DLA 
is  beyond  the  scope  of  this  paper.  Let  us  however  emphasize  that  although 
DLA  do  not  have  exact  affine  self-similarity  properties,  our  maxima  voting 
algorithm  recovers  enough  statistical  information  about  the  renormalization 
parameters,  to  test  physical  models  of  these  aggregates.  We  believe  that 
such  robust  voting  schemes  can  be  used  to  test  a  large  class  of  renormaliza¬ 
tion  models  for  physical  multifractals. 
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Figure  8:  (a)  Center  region  of  a  DLA  cluster  of  108  particles,  which  is  in¬ 
cluded  in  a  circle  of  fixed  radius,  (b)  measure  showing  the  intersection  points 
of  a  DLA  cluster  and  a  circle  of  fixed  radius  (data  provided  by  Arneodo). 


Figure  9:  This  curve  gives  the  decay  of  log  N(s)  as  a  function  of  log. s,  where 
JV'(s)  is  the  number  of  wavelet  transform  modulus  maxima  at  the  scale  s  for 
50  DLA  measures. 
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Figure  10:  (a)  Distribution  of  votes  in  the  plane  (log  p ,  log  /)  of  the  parameter 
space,  for  50  DLA  measures.  The  votes  are  distributed  along  a  straight  line 
of  equation  logp  =  0.63  log  I.  (b)  Distribution  computed  for  50  stochastic 
Fibonacci  Cantors. 
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Appendix:  proof  of  Theorem  1 


The  wavelet  transform  given  in  equation  (5)  can  be  rewritten  a  convolution 
product 

Wf(s,x)  =  f*-tps(x),  (21) 


where  ips(x)  =  -).  The  wavelet  transform  of  a  fractional  Brownian 

motion  is  a  stationary  Gaussian  process  [8]  with  a  power  spectrum  P(o>) 
equal  to 


PH 


<72|-0(,H|2 
M2ff+1  5 


(22) 


where  ^(u>)  is  the  Fourier  transform  of  ip(x).  The  density  of  m  of  modulus 
maxima  of  a  function  is  related  to  the  density  e  of  local  extrema  and  its 
density  2'  of  zero-crossings  by 


e  +  z 


The  density  of  zero-crossings  of  a  differentiable  Gaussian  process  whose  au¬ 
tocorrelation  is  R(t )  is  [13] 


2  = 


(24) 


where  p(”)(r)  is  the  nth  derivative  of  R(t).  The  density  of  extrema  of  a 
differentiable  process  is  equal  to  the  density  of  zero-crossings  of  the  derivative 
of  the  process.  If  the  autocorrelation  of  a  process  is  P(  r),  the  autocorrelation 
of  its  derivative  is  —  R(2\t).  We  thus  derive  from  equation  (24)  that  for  a 
stationary  Gaussian  process  that  is  twice  continuously  differentiable,  the 
density  of  extrema  is 

p (25) 

V  -IT*  JJM(O)'  ( 


The  derivatives  of  the  autocorrelation  function  in  zero  is  related  to  the  power 
spectrum  by 


«0)  =  j 

r+oo 

j  P(io)dto, 

—  oo 

(26) 

-  Pt2)(0)  = 

f- 1-0© 

/  to2  P(u;)duj, 

J  —  OO 

(27) 

and 

P(4)( 0)  =  j 

f  +  OO 

1  LO4P(u))di0. 

—  OO 

(28) 
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If  the  wavelet  is  twice  continuously  differentiable,  the  convolution  of  a  Frac¬ 
tional  Brownian  motion  with  a  dilated  wavelet  is  a  stationary  Gaussian  pro¬ 
cess  that  is  twice  differentiable.  Inserting  equation  (22)  in  equations  (26-28) 
yields 

/-{-oo 

M-1-2Hj#(")|2  dLV,  (29) 

-co 

/+CO 

|u;|1-2i?|#(w)|2  duj,  (30) 

-co 

/+CO 

\lu\3-2H\4>(lo)\2  du>.  (31) 

-co 

As  a  consequence  of  equation  (23),  the  average  density  of  modulus  maxima 
of  a  wavelet  transform  is  given  by  the  average  density  of  its  zero-crossings 
and  of  its  extrema.  From  equations  (24-25)  and  (29-31),  we  derive  that  the 
average  density  of  modulus  maxima  of  the  wavelet  transform  at  the  scale  s 
of  a  fractional  Brownian  motion  is 


with 


and 
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This  finishes  the  proof  of  Theorem  1. 
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